Resistivity of inhomogeneous quantum wires 
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We study the effect of efectron-electron interactions on the transport in an inhomogeneous quan- 
tum wire. We show that contrary to the weil-known Luttinger iiquid resuit, non-uniform interactions 
contribute substantially to the resistance of the wire. In the regime of weakly interacting electrons 
and moderately low temperatures we find a linear in T resistivity induced by the interactions. 
We then use the bosonization technique to generalize this result to the case of arbitrarily strong 
interactions. 

PACS numbers: 71.10.Pm 



Since the first measurements of quantized dc conduc- 
tance in quantum wires [H , transport properties of these 
systems have generated a lot of interest. Theoretically, 
such one-dimensional conductors cannot be described by 
the conventional Fermi liquid theory, but rather form a 
ualitatively different state known as the Luttinger liquid 
Recently several characteristic signatures of the Lut- 
tinger liquid state have been reported experimentally in 
quantum wires 0, 0| • Perhaps even more interestingly, a 
number of experiments show anomalies in the transport 
properties of these systems in the form of a small struc- 
ture below the first plateau of quantized conductance 
[1, @, 0, @, S [loll . which are not expected in the Luttinger 
liquid theory The so-called "0.7-structure," which 

develops at finite temperature in low-density wires, is 
an example of such deviations from perfect quantization 
BUSES]- While its precise origin still remains unclear, 
this feature is most likely related to correlations between 
electrons, which initiated various attempts to study the 
effect of interactions on the electronic transport in these 
devices [II [IS [II]. 

A number of recent theory papers [H , 14 1 studied the 
model of a quantum wire device in which interactions are 
present only in a small region of a one-dimensional elec- 
tron system between two non- interacting leads. If the 
size of the interacting region does not significantly ex- 
ceed the Fermi wavelength of the electrons in the wire, 
the interactions give rise to backscattering of either sin- 
gle electrons or pairs, resultin g in significant corrections 
to the quantized conductance [l4j . On the other hand, if 
the interaction strength varies smoothly over a long dis- 
tance, such backscattering processes are expected to be 
exponentially weak and can be neglected. In this regime 
a model of non-uniform Luttinger liquid, with parame- 
ters gradually varying as a function of position is appro- 
priate. Studies of such a model found no correction to 
the quantized dc conductance of the wire [ill ] . It is thus 
natural to conclude that inhomogeneities of interacting 
quantum wires at large scales d 3> kp 1 do not affect the 
dc transport beyond the exponentially small backscatter- 
ing corrections. 

In this paper we show that even at kpd ^> 1, when 



the backscattering processes [14j can be ignored, the in- 
homogeneity of the interaction strength in the wire gives 
rise to a finite resistivity at non-zero temperature. 

We start by considering an infinite one-dimensional 
system of weakly interacting spinless electrons, with 
quadratic dispersion e p = p 2 /2m. In this simple model, 
the electron density n is assumed to be uniform, but the 
strength of the electron-electron interactions varies along 
the wire. We describe these inhomogeneous interactions 
by the potential 



V{x,y) = V(x -y) rj 



y 



(i) 



Here V(x — y) is the conventional electron-electron re- 
pulsive interaction. Coulombic in nature, it is screened 
by the nearby gates, and for simplicity we will treat it 
as a short-range interaction. The non-uniformity of the 
system is then encoded in the dimensionless function r], 
which varies at a length scale d, large compared with 
both the Fermi wavelength and the range of the interac- 
tion potential V(x — y). 

In order to compute the resistance of the wire, we en- 
force a dc current / to flow through the system. The 
electrons in the wire then acquire a drift velocity propor- 
tional to this current: Vd — I /ne. In the reference frame 
moving along the wire with velocity Vd the electronic sub- 
system is in equilibrium, as pointed out by Pustilnik et al. 
[13 ] in the context of Coulomb drag between two paral- 
lel wires. This equilibrium is characterized by a Fermi 
energy and a temperature T. 

When viewed in the stationary reference frame, where 
the electric current does not vanish, the electrons are 
no longer in thermodynamic equilibrium. In particu- 
lar, their occupation probabilities cannot, in general, be 
expressed as a Fermi function of the energy. However, 
at T <C €f the occupation probabilities of the left- and 
right-moving states near the Fermi level can still be ap- 
proximated by Fermi functions, albeit with two different 
temperatures, Tl and Tr. To show that, we note that 
the electron energy e p changes to a different value e p in 
the stationary frame. Considering a state p near the right 
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FIG. 1: Electronic dispersion in the stationary frame, with 
an example of non momentum-conserving process. The right- 
and left- movers Fermi energies read e^' L = ef(1 ± 2vd/vp). 



Fermi point pp, to first order in p — pp, we have 



e p - e F _ vp(p-p F ) 



l p -lp vp(p~p F ) 



vf 



where in the stationary frame vp = vp + Vd, p = p + mvd, 
Pp — pp + mvd- Thus, the occupation probability of this 
state can be expressed in terms of its energy e p in the 
stationary frame as 

1 1 



e ( £p -£ F )/T + i e (i P -a F )/T R + i ' 

where Tr — Tvf/vf = T(l + v c i/vp). Similarly, for the 
electrons near the left Fermi point — pp one finds the 
occupation probability given by the Fermi function with 
the effective temperature Tp = T(l — Vd/vp). 

One may expect that the electron-electron interactions 
induce thermalization between right- and left-movers 
through two-particle scattering processes like the one 
shown in Fig. [T] In a uniform system, two-particle pro- 
cesses cannot lead to thermalization [16fl because of the 
conservation of both energy and momentum: the interac- 
tions either exchange the momenta of the two electrons 
or leave them unchanged. However, in our model the 
non-uniformity of the interaction potential (Q]) breaks the 
translational invariance of the system, and allows for two- 
particle scattering processes that conserve energy but not 
momentum. 

A typical two-particle process shown in Fig. [1] de- 
scribes the scattering from a state with momenta (p, k) to 
(p', k'), and is accompanied by an overall loss of momen- 
tum. Since this process involves a transfer of energy from 
the "warmer" right-moving branch to the "colder" left- 
moving one, it is expected to occur more frequently than 
the inverse process (p',k') — > (p, fc), so that on average, 
the two subsystems lose more momentum than they gain. 
Note that unlike Ref. Q, the typical change of momen- 
tum for the processes shown in Fig. [T] is small compared 
to the Fermi momentum, and the rate of such processes 
will not become exponentially small at kpd^> 1. 

The decrease in momentum can be interpreted as a 
result of a damping force acting on the electrons. To 



maintain constant current, it has to be balanced by a 
driving force, which stems from a local electric field, gen- 
erated as a response of the system to the current bias [15| ■ 
Since the damping force is proportional to the tempera- 
ture difference between the two subsystems, Tr — Tl cx /, 
this local electric field is proportional to the applied cur- 
rent, and the proportionality coefficient is defined as the 
resistivity. 

Let us compute the resistivity at temperatures T 3> 
hvp/d. We isolate a small segment of wire taken at posi- 
tion x, with the length Ax in the range hvp/T <C Ax -C 
d. The driving force eE(x)nAx acting on this segment 
of wire as a result of the local electric field E(x) = p(x)I 
compensates for the damping force AF due to the inter- 
actions, so that the resistivity can be written as 



p(x) 



AF 

enlAx 



(2) 



We compute the damping force as the change in momen- 
tum per unit time, using the Fermi golden rule 



AF = — \V pk -p'k'\ 2 S(e p + e k - e p > - e k >) 

x(p' + k'-p-k) f*tf(l ~ /«)(! - /£).(3) 



Here V p k- P 'k' is the matrix element of the interaction po- 
tential for scattering from an initial state (p, k) to a final 
state (p\ k'), as shown in Fig. [1] the occupation numbers 
f ' L are given by the Fermi-Dirac distribution evaluated 
with the appropriate temperatures Tr ; j,. One can easily 
check that at Tr = Tp expression ([3]) vanishes. Then in 
the linear response regime, one can expand the occupa- 
tion numbers f R and f L to first order in Tr — Tp oc /, 
and find that A_F is proportional to the applied current. 

To first order in the interaction potential, the matrix 
element V p k-, P 'k' is given by: 

W*' = [ X+AX dye* P ' +k '* P ~ k y [1/(0) - V(2k F )} V (y), 

(4) 

where V(0) and V(2kp) are respectively the zero and 
2kF Fourier components of the interaction potential in- 
troduced in Eq. (l}. Substituting Eq. g} into Eq. (J3j) , 
one readily sees that a constant value of r\ enforces the 
conservation of momentum and leads to a vanishing re- 
sult. The dominant non- vanishing part thus involves the 
gradient of rj, and contributes to AF as (d x r)) 2 . 

Performing the remaining momentum summations, the 
resistivity evaluated to second order in the interaction 
then takes the form 



p(x) 



h T 
64e 2 neF 



V(0) - V{2k F )\ 2 fdr 1 {x)'" 2 



irhvF 



\ dx 



(5) 



This result was obtained in the regime of temperatures 
T 3> hvp/d, for which the position integrals coming from 
the matrix element Q could be easily simplified. 



3 



Our method provides a clear physical picture of the 
origin of the resistivity, in the simple case of weakly in- 
teracting spinless fermions. We now turn to the case 
of arbitrarily strong interactions, where we derive the 
expression for the resistivity using a bosonized Hamil- 
tonian. In addition, we account for electron spins and 
allow for a non- uniform density n(x), as a result of the 
surrounding gates and impurities in the substrate. 

Following Ref. [II], we generalize the Tomonaga- 
Luttinger model of interacting one-dimensional electron 
systems to account for inhomogeneities by allowing for 
position dependence of the Luttinger-liquid parameters. 
In the case of electrons with spins, this procedure yields 



H 



Hp + H a 



tlx 



dx 



hu p {x) 

%u a (x) 
2tt 



K p {x) {d x 6 P Y + 



K a {x) (d x 9 a ) 



■2 , {dx4> P Y 



K p {x) ^ 

(d x <j>gf 



(6a) 
(6b) 



2 9<J (x) f 
ax k cos 



[27ra(x)Y 



[2^2^ 



(6c) 



Here the short distance cutoff a is assumed to be a 
function of x. In the limit of a homogeneous system, 
kpd — > oo, the coupling constant g a renormalizes to zero 
at large length scales; at the same time, the parameter 
K a approaches unity as K a = 1 + g cr /2ir?m cr , as required 
by the SU(2) symmetry Q. We assume that kpd is suf- 
ficiently large for the system to be near this limit at ev- 
ery point x, i.e., g rT (x)/Uu rT (x) <C 1. Previous works on 
the inhomogeneous Luttinger liquid model were either 
restricted to spinless electrons [ljj, or discarded X7| the 
cosine term in Eq. (|6c[) . invoking the irrelevance of g a at 
low energies. In our case its contribution to the resistiv- 
ity is as important as that of the quadratic part of H a . In 
the derivation below, we assume that the inhomogeneities 
of the system are weak, e.g. \u p (x) — u p (0)\ <§; u p (0). 

The resistivity can now be computed following a 
method similar to the one outlined in Ref. [l3| in the con- 
text of a quantum wire in the Wigner crystal regime. As 
one applies an electric current I = Iq cosut, the electrons 
start moving in the wire. In the dc limit ui — > 0, we can 
assume that all electrons move in phase, so that at time 
t their position has shifted by a distance proportional 
to the injected charge q(t) = J w _1 smut. As a conse- 
quence, we need to evaluate all the position-dependent 
parameters in the Hamiltonian ([6]) at the true time- 
dependent position of the electrons, which amounts to 
replacing x — > x + q(t)/en(x). In the regime of linear re- 
sponse, we only need to expand the Hamiltonian to first 
order in q(t), 



H = dx[H( 



en{x) 



H'(x) 



(7) 



where 7i(x) is the Hamiltonian density when no cur- 
rent is applied, and H'(x) is obtained from Ti(x) by re- 
placing the position-dependent parameters g a {x) / a(x) 2 , 
Uv{x)K v (x) and u v (x) / 'K v {x) (where v — p, a) by their 
derivatives with respect to x. 

In the conventional Luttinger-liquid theory, the cur- 
rent I = q is usually viewed as an excitation of the charge 
mode, and q{t) thus appears as a dynamical variable pro- 
portional to (j) p . Then the linear in q part of Eq. cor- 
responds to cubic terms such as (f>p(d x (f>p) 2 ■ These cubic 
terms are usually disregarded as irrelevant perturbations 
to the Luttinger liquid Hamiltonian. Nevertheless, the ef- 
fect of such perturbations should be addressed, because 
without them no contribution to transport arises from a 
non-uniform interaction [H| . In what follows, it will be 
more convenient to treat q{t) as an external parameter. 

The oscillatory perturbation in the Hamiltonian ([7]) 
acts as an external driving force, which leads to the cre- 
ation of spin and charge excitations and dissipation of 
the energy from the driving force to the wire. The en- 
ergy W dissipated into these excitations in unit time may 
be obtained using the Fermi golden rule. In the limit of 
weak applied current, it is expected to be quadratic in 
the amplitude Iq of the current oscillations. This allows 
us, by comparison with the Joule heat law W = IqR/2, 
to deduce the expression for the resistance R of the wire. 
Then in the dc limit w->0wc find 



R = 



e 2 h 



dx dy 
n(x)n(y) 



dtt{H'(x,t)H'(y,0)), (8) 



where (. . .) corresponds to the thermodynamic average. 
The last integral in Eq. ((5J) falls off rapidly at \x — y\ ^S> 
hu Pt(T /T. Thus at T ^> hu Pt(J /d, we can reduce the ex- 
pression ([5]) to a single integral in space, whose integrand 
we identify with the resistivity p(x) of the wire. 

As both the charge and spin modes dissipate energy 
throughout the wire, the total resistivity is given by the 
sum of a charge and a spin contribution, p(x) — p P (x) + 
Pa{x), which can be computed separately. Substituting 
the charge Hamiltonian (|6b[) in Eq. ([5]), and performing 
the remaining time integral, one can extract the charge 
contribution to the resistivity 



P P ( X ) 



T 



d x K p {x) 
lire 2 %Up(x)[n(x)] 2 \ K p (x) 



(9) 



This result holds at temperatures in the range hu p /d <C 
T <C Dp, where the charge bandwidth D p ~ Tinu p . 

One can use the expression @ to recover our earlier 
result © for weakly interacting spinless electrons. In 
this case, upon bosonization the Hamiltonian of the sys- 
tem takes a form equivalent to H p with u p (x) —> vf and 
Kp(x) -> 1 - [V(0) - V(2k F )] i](x)/2Trhv F . Substituting 
these expressions into Eq. |9]) and expanding to second 
order in the interaction, one reproduces the result ([5|). 

The spin contribution to the resistivity consists of two 
terms arising from substituting the quadratic part and 
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the cosine part of the spin Hamiltonian in Eq. {SJ. The 
contribution of the quadratic part of H a can be ob- 
tained from Eq. @ by replacing the charge parame- 
ters with their spin counterparts. This result is further 
simplified by using the low-energy expansion K a (x) — 
1 + y a {x)/2, where we introduced the dimensionless pa- 
rameter y a {x) = g ( y{x)/i:fiu cr {x). 

In the case of weakly interacting electrons, the correc- 
tion to the parameter K a in the quadratic part of Eq. (|6c|) 
accounts for the gi-coupling of the z components of the 
electron spins, while the cosine term in H a is associated 
with the remaining x and y components Q. Then from 
the spin symmetry of the system, the cosine term of the 
spin Hamiltonian should contribute twice as much to the 
resistivity as the quadratic part. We expect this result 
to hold for arbitrarily strong interactions; in the case of 
weakly interacting electrons, this conclusion can easily 
be verified [l8[. Combining these two terms, the spin 
contribution to the resistivity reads: 



Pa{x) 



lib 



T 



32ne 2 hu a {x)[n{x)] 



:[d xVa {x)] 2 . (10) 



Again, we restricted ourselves to the range of moderately 
low temperatures, Tiu a /d <C T <C D a , where the spin 
bandwidth is given by D a ~ Tinu a . 

The comparison of the two contributions © and (jTUJ) 
to the resistivity of the wire suggests the strongest effect 
in the regime of low electron density, when the electron 
correlations are strong. In this case the exchange cou- 
pling J of electron spins, which sets the spin bandwidth 
D(j , is strongly suppressed, so that D a <C D p . As a re- 
sult, we expect the spin part (flOj) of the resistivity to 
be the dominant contribution in this regime, due to the 
reduced spin velocity. It is worth pointing out that y a 
is only marginally irrelevant, so that while it renormal- 
izes to zero at low temperature, it does so logarithmi- 
cally, as 1/ \og(D a /T). Furthermore, most experimental 
measurements are carried out at fixed temperature, while 
varying the electron density in the wire. In this config- 
uration, the logarithmic dependence of y a suggests that 
this parameter increases as the interaction in the wire 
becomes stronger. 

Our results are relevant to experiments on wires longer 
than the length l eq associated with the processes of equi- 
libration in the moving frame. In shorter wires, with 
length L <C l e q, we expect the resistivity to be sup- 
pressed by an additional factor of order L/l eq . This 
raises a fundamental question of the equilibration in a 
one-dimensional system of interacting electrons. In the 
weakly interacting case, it is believed that the leading 
equilibration mechanism is due to three-particle collisions 
[161 ] involving states near the bottom of the electronic 
band. One expects such processes to be strongly sup- 
pressed at low temperatures, corresponding to a large 



l eq . However, such a treatment is not applicable beyond 
the limit of weak interactions. While we expect stronger 
interactions to make thermalization easier, a detailed in- 
vestigation of the equilibration processes will be neces- 
sary to access the full temperature dependence of the 
resistivity in this regime. 

In summary, we have shown that the interactions be- 
tween electrons in a long inhomogeneous quantum wire 
give rise to a finite resistivity p — p p + p a , given by 
Eqs. ([9]) and (fTQ| . This resistivity is due to the weak vio- 
lation of the momentum conservation in electron-electron 
collisions, caused by the inhomogeneities on long spatial 
scales d^> kp . Our results can be tested experimentally 
by measuring the temperature and density dependences 
of the resistance of long quantum wires. 
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